An Unconditionally Energy Stable Penalty Immersed Boundary Method for Simulating the Dynamics of an Inextensible Interface Interacting with a Solid Particle
نویسندگان
چکیده
In this paper, a novel penaltymethod based on the immersed boundary formulation is proposed for simulating the transient Stokes flow with an inextensible interface enclosing a suspended solid particle. The main idea of this approach relies on the penalty techniques by modifying the constitutive equation of Stokes flow toweaken the incompressibility condition, relating the surface divergence to the elastic tension σ to relax the interface’s inextensibility, and connecting the particle surface-velocity with the particle surface force F to regularize the particle’s rigid motion. The advantage of these regularized governing equations is that when they are discretized by the standard centered difference scheme on a staggered grid, the resulting linear system can easily be reduced by eliminating the unknowns ph, σh and Fh directly, so that we just need to solve a smaller linear system of the velocity approximation uh . This advantage is preserved and even enhanced when such approach is applied to the Po-Wen Hsieh was partially supported by the National Science Council of Taiwan under the Grant NSC 102-2115-M-033-007-MY2. Ming-Chih Lai was partially supported by the National Science Council of Taiwan under the Grant NSC 101-2115-M-009-014-MY3. Suh-Yuh Yang was partially supported by the National Science Council of Taiwan under the Grant NSC 101-2115-M-008-008-MY2. P.-W. Hsieh Department of Applied Mathematics, Chung Yuan Christian University, Jhongli City 32023, Taoyuan County, Taiwan e-mail: [email protected] M.-C. Lai Center of Mathematical Modeling and Scientific Computing, Department of Applied Mathematics, National Chiao Tung University, Hsinchu 30010, Taiwan e-mail: [email protected] S.-Y. Yang (B) · C.-S. You Department of Mathematics, National Central University, Jhongli City 32001, Taoyuan County, Taiwan e-mail: [email protected] C.-S. You e-mail: [email protected] 123 Author's personal copy
منابع مشابه
Simulating the dynamics of inextensible vesicles by the penalty immersed boundary method
In this paper, we develop an immersed boundary (IB) method to simulate the dynamics of inextensible vesicles interacting with an incompressible fluid. In order to take into account the inextensibility constraint of the vesicle, the penalty immersed boundary (pIB) method is used to virtually decouple the fluid and vesicle dynamics. As numerical tests of our current pIB method, the dynamics of si...
متن کاملAn immersed boundary method for simulating vesicle dynamics in three dimensions
We extend our previous immersed boundary (IB) method for 3D axisymmetric inextensible vesicle in Navier-Stokes flows, [W.-F. Hu, Y. Kim, M.-C. Lai, An immersed boundary method for simulating the dynamics of three-dimensional axisymmetric vesicles in Navier-Stokes flows, J. Comput. Phys. 257 (2014)] to general three dimensions. Despite a similar spirit in numerical algorithms to the axisymmetric...
متن کاملA Fractional Step Immersed Boundary Method for Stokes Flow with an Inextensible Interface Enclosing a Solid Particle
In this paper, we develop a fractional step method based on the immersed boundary (IB) formulation for Stokes flow with an inextensible (incompressible) interface enclosing a solid particle. In addition to solving for the fluid variables such as the velocity and pressure, the present problem involves finding an extra unknown elastic tension such that the surface divergence of the velocity is ze...
متن کاملUnconditionally Energy Stable Immersed Boundary Method with Application to Vesicle Dynamics
We develop an unconditionally energy stable immersed boundary method, and apply it to simulate 2D vesicle dynamics. We adopt a semi-implicit boundary forcing approach, where the stretching factor used in the forcing term can be computed from the derived evolutional equation. By using the projection method to solve the fluid equations, the pressure is decoupled and we have a symmetric positive d...
متن کاملDynamics of Love-Type Waves in Orthotropic Layer Under the Influence of Heterogeneity and Corrugation
The present problem deals with the propagation of Love-type surface waves in a bedded structure comprises of an inhomogeneous orthotropic layer and an elastic half-space. The upper boundary and the interface between two media are considered to be corrugated. An analytical method (separation of variables) is adapted to solve the second order PDEs, which governs the equations of motion. Equations...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
عنوان ژورنال:
- J. Sci. Comput.
دوره 64 شماره
صفحات -
تاریخ انتشار 2015